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Analysis  of  Scattering  and  Radiation  of  Acoustic 
Waves  From  an  Axisymmetric  Structure  in  an 
Infinite  Fluid  Medium 


Introduction 

Determination  of  the  acoustic  radiation  from  a  vibrating  surface  in  an  infinite,  inviscid  fluid 
medium  requires  solution  both  of  the  equations  of  motion  of  the  structure  and  of  the  Helmholtz 
equation  for  the  velocity  boundary  condition  normal  to  the  structure’s  surface.  However,  the 
availability  of  an  exact  solution  of  the  Helmholtz  boundary  value  is  usually  precluded  because  most 
structures  of  interest  have  a  complicated  surface  geometry.  As  a  result,  the  solution  of  acoustic 
radiation  problems  has  generally  employed  two  alternative  approaches. 

In  the  first  approach1  the  structure  is  contained  within  a  sphere,  and  everything  within  the  sphere  is 
regarded  as  the  structure.  The  surrounding  fluid  is  idealized  by  fluid  finite  elements,  and  the  structure 
is  idealized  by  structural  finite  elements.  The  exact  solution  of  the  Helmholtz  equation  for  a  spherical 
cavity  in  an  infinite  medium  is  available,  and  this  solution  is  coupled  with  the  solution  for  the  fluid- 
structure  finite  element  sphere. 

In  the  second  approach2  the  integral  equation  solution  of  the  Helmholtz  boundary  value  problem 
is  derived.  This  “surface  Helmholtz  integral” is  derived  from  the  conventional  Helmholtz  integral  for 
a  field  point  exterior  to  the  cavity  by  a  limiting  process  in  which  the  field  point  approaches  the  sur¬ 
face. 

Use  of  the  first  approach  requires  solution  of  a  large  fluid-structure  finite  element  problem 
whenever  a  long,  thin  structure  is  to  be  analyzed.  Use  of  the  second  approach  yields  a  nonunique 
solution  in  the  case  where  the  excitation  frequency  is  within  a  certain  bandwidth  around  the  natural 
frequency  of  the  interior  Helmholtz  problem.  Two  distinct  methods  3  4  are  available  to  alleviate  these 
difficulties.  Although  they  are  not  detailed  here,  these  methods  are  incorporated  in  the  computer 
program. 

In  earlier  works2-3-5  on  the  solution  of  the  Helmholtz  surface  integral,  it  is  assumed  that  pressure 
and  velocity  over  each  “finite  surface  element”  are  constant;  that  is,  each  finite  surface  element  acts 
as  a  piston.  Our  analysis  gives  the  relation  between  pressure  and  velocity  at  the  centroid  of  each  finite 
surface  element.  The  total  set  of  finite  surface  elements  constitutes  the  total  fluid  surface  area  in 
contact  with  the  structure. 

In  the  earlier  works,  the  response  of  the  structure  (normally  given  at  the  nodes)  and  the  pressure- 
velocity  relation  obtained  at  the  centroids  of  the  finite  surface  elements  are  made  compatible  by  in¬ 
terpolating  the  response  at  the  node  points  from  the  response  at  the  centroids,  However,  since  most 
structures  are  analyzed  using  the  finite  element  method,  it  seems  appropriate  to  divide  the  fluid 
surface  and  the  structural  surface  in  contact  with  each  other  by  the  same  network  of  grid  lines.  In  this 
way,  each  of  the  finite  surface  elements  has  a  counterpart  structural  finite  element  on  the  structural 
surface.  (Previous  works  do  not  impose  such  a  restriction.)  We  can  then  further  require  that  the 
velocity  and  pressure  distribution  on  both  sets  of  elements  be  consistent.  In  essence,  then,  we  choose 
to  view  the  fluid-structure  interaction  problem  as  simply  an  extension  of  the  finite  element  analysis  of 
structures.  This  is  the  new  feature  examined  in  this  study. 
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The  results  obtained  thus  far  show  that,  for  a  given  criterion  of  accuracy,  the  number  of  finite 
surface  elements  required  to  solve  the  Helmholtz  surface  integral  is  less  than  that  required  for 
structural  analysis.  This  assures  us  that  the  finite  element  model  of  the  structure  and  its  counterpart  in 
fluid  is  sufficient  for  determining  structural  response  in  an  infinite  fluid  medium. 

Although  the  algebraic  development  of  the  work  here  is  restricted  to  an  axisymmetric  case,  the 
method  can  easily  be  extended  to  the  nonaxisymmetric  case. 

Analysis 

The  following  analysis  is  divided  into  three  parts.  The  first  part  deals  with  the  numerical  solution  of 
the  Helmholtz  surface  integral,  the  second  describes  the  equations  of  motion  of  the  structure  in  the 
fluid,  and  the  third  deals  with  the  response  of  the  structure  and  scattering  of  the  incident  wave  train 
along  the  axis  of  symmetry. 

Solution  of  Helmholtz  Surface  Integral 

Radiation  from  an  inclusion  in  an  infinite  fluid  medium  is  described  by  the  boundary  value 
problem  where  the  Helmholtz  equation  represents  the  field  equation: 

(V2+k2 )p(x)  =0  (1) 

subject  to  the  Sommerfeld  radiation  condition  at  infinity: 

*<!s  -  “■»  ■  °  ® 

and  the  pressure  p(x ')  and  velocity  vn(x.')  normal  to  the  surface  prescribed  on  the  inclusion  surface. 

The  solution  of  the  above  boundary  value  problem  is  given  by  the  Helmholtz  integral 

ikR  iRR 

=  h  f  {p(^  h  (5Lir)  -  ikp*c,\  <*>  >  dA,  o) 

A 

where  (see  figure  1): 

x*.  i define  the  vector  from  the  origin  of  the  coordinate  system  to  the  field  point  A  and  the  vector  in 
the  fluid  to  the  general  surface  point  B,  respectively 

n  is  the  exterior  normal  to  the  surface  at  point  B 


R  =  |j*-a| 


V(x)  =  velocity  vector  of  point  B 

V„(x)  =  V(x)*fl  =  component  of  V(x)  normal  to  the  surface 


k  =  o>/c 

V1  =  Laplacian  operator 

2 


TR  5821 


u>  =  frequency  of  excitation  in  radians/sec 
c  =  velocity  of  sound  in  the  fluid 
e  =  mass  density  of  the  fluid . 


A 


Figure  1.  Schematic  for  Helmholtz  Integrals 


The  Helmholtz  integral  (equation  (3))  gives  the  pressure  at  a  field  point  A  for  the  case  in  which  the 
velocity  and  pressure  distribution  on  the  inclusion  surface  are  given.  The  problem  addressed  here  is— 
given  a  velocity  distribution  on  the  inclusion  surface,  find  the  corresponding  pressure  distribution  on 
the  surface.  This  was  obtained  from  equation  (3)  by  considering  the  limiting  case  where  the  field  point 
A(x ' ')  approaches  the  surface  point  C(x ').  This  case  is  considered  in  detail  by  Chertock2,  where  the 
resulting  expression,  called  the  Helmholtz  surface  integral,  is  an  integral  equation  for  p(x/)  when 
Vn(x)  is  prescribed: 

f*  ~  ikR  ikR 

P(XI)  =  ^  j  (p(x)  £  — R-  }  dA  (4) 


where  the  asterisk  (*)  over  the  integral  means  that  the  principal  value  of  the  improper  integral  is  to  be 
taken  on  the  total  surface  A.  Note  that  as  R-’-O,  the  integrand  exp(ikR)/R-*°°,  which  accounts  for  the 
difference  between  the  integrals  of  equations  (3)  and  (4). 

The  actual  shapes  of  structures  normally  encountered  rarely  conform  to  the  shapes  (sphere, 
cylinder,  etc.)  for  which  the  exact  solution  of  equation  (4)  is  available.  As  a  result,  numerical  analysis 
in  which  the  surface  is  divided  into  a  set  of  finite  surface  elements  is  used.  Assume  that  the  total  fluid 
surface  is  divided  into  J  number  of  “finite  surface  elements”;  then  the  pressure  at  a  grid  point  C(x ')  is 
given  by 

P(?‘>  “  h  2  /  {p%)  Tn  (^1")  "  ikpcVn(V  1  dAi-  (5> 

j*l  Aj 
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If  the  pressure  and  velocity  over  each  “finite  surface  element”  are  regarded  as  constant  (which  is  the 
usual  assumption),  p(xj)  and  Vn0cp  will  drop  out  of  the  integral  sign,  so  that 


1  f  3  ikR  f  ikR 

'l*l)  =  2v  2  {P^j>  J  9n  >  dAj  -  Vn(^j)j  ikpr  r  dA’: 

J=1  Aj  J  Aj 


(6) 


Note  that  \  '  is  the  centroid  of  a  particular  finite  surface  element  and  R  =  |x  '-*j|.  Also,  it  should  be 
clear  that  both  of  the  integrals  in  equation  (6)  will  have  to  be  evaluated  numerically. 


As  mentioned  earlier,  the  method  under  consideration  here  employs  the  consistent  distribution 
approach.  In  the  finite  element  method,  the  displacement  at  a  point  within  a  particular  plate  finite 
element  is  described  by  a  basis- function  (normally  a  polynomial  in  local  coordinates  of  the  order  such 
that  the  number  of  coefficients  of  the  polynomial  equal  the  total  number  of  variables  available  at  the 
total  number  of  grid  points  describing  the  finite  element).  Consider,  for  example,  the  triangle 
described  by  three  grid  points  1, 2, 3  (see  figure  2). 


Y 


The  terms  U,,  U4,  U7  (the  displacements  normal  to  the  plane)  and  the  two  rotations  at  each  grid  point 
that  can  produce  the  normal  displacement  of  the  plate  account  for  the  total  of  nine  variables. 
Utilization  of  these  nine  variables  allows  us  to  select  a  polynomial  in  two  dimensions  (x,y — which  is 
the  local  coordinate  system  of  the  triangular  plate)  with  nine  coefficients: 

uz  =  Cj  +  C2X  +  <:3y  +  (:4X  +  <-5xv  +  C  \  +  CjX '  +  ( x'  v+xv  )  +  C  v  '  . 

These  nine  coefficients  can  be  expressed  in  terms  of  .nine  known  quantities.  Without  deriving  the 
expressions  explicitly  for  the  coefficient  Cj,  we  will  write  uz  in  matrix  notation: 


4 
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u 

7 


.1 .  ( x ,  v )  u  . 


i  =  1 


where  a^x.y)  are  the  polynomials  in  local  coordinates  x,y,  and 


(7) 


Selection  of  the  basis  function  for  an  axisymmetric  conical  shell  element  is  detailed  in  appendix  A. 


Note  that  the  nonlinear  variation  of  the  displacement  in  the  previous  example  requires  not  only  the 
displacements  at  the  grid  points,  but  also  the  rate  of  change  of  the  displacements.  Since  the  general 
purpose  computer  programs  based  on  the  finite  element  method  always  give  the  translational  and 
angular  velocities  at  the  grid  points,  the  use  of  nonlinear  distribution  does  not  require  any  change  to 
the  existing  computer  programs.  This  is  not  the  case  with  the  pressure  distribution  problem.  Equation 
(5)  gives  pressure  information  only,  whereas  for  the  nonlinear  pressure  distribution  case,  we  also  need 
the  rate  of  change  of  pressure  at  the  grid  points.  To  obtain  this  information,  we  will  differentiate 
equation  (5)  with  respect  to  the  local  spatial  variables: 


3 


3 1 


k 


p(x  !  ) 


* 


where 


ikpcV  (x  . )  1  1  it  A  j  ,  (8) 

n  ]  ot:  K 

k 


ikR 


k  =  1,2 

n, ,  1 1 ,  t2  form  a  local  right-hand  coordinate  system  at  point  C(x ') 
n  is  the  normal  to  the  surface  at  C(x ') 

tj,  t2  are  tangent  to  the  surface  at  C(x')  and  are  at  right  angles  to  each  other.  Either  t,  or  t2 
(whichever  is  convenient)  is  fixed  to  define  the  local  coordinate  system  (see  x,y,z  coordinate  system  in 
figure  2);  the  remaining  axis  should  form  the  right-hand  system. 


Then,  with  both  p(x ')  and  3/(3tp)*p(x ')  known  at  the  grid  points,  an  interpolation  function 
(corresponding  to  equation  (7))  can  be  written  for  the  pressure. 


In  accordance  with  the  notations  used  in  equation  (7)  for  the  consistent  distribution,  let  us  define 
pressure  and  velocity  for  the  jth  “finite  surface  element”: 


T  & 

p  ( x  .  )  =  .1  P 
"i  i  i 


T  5 

V  (X.)  =  V 
n  -)  -)  -j  > 


(9) 
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where  I  refers  to  the  variable  defined  in  local  coordinates. 

In  the  analysis  we  also  considered  the  second  case  in  order  to  see  how  much  difference  (if  any)  the 
linear  distribution  of  pressure  over  the  “finite  surface  elements”  would  make.  In  the  linear  case,  the 
first  set  of  equations  of  equation  (9)  will  change : 

T  ? 

p(x  .  )  =  c  .  P 
1  .1  1 

V  (x.)  =  aT  VZ  (10) 

n  ".I  J  M 

To  increase  notational  flexibility,  we  will  use  equation  (10)  for  both  cases,  with  the  understanding  that 

for  the  consistent  distribution  case  qJ  =  qJ  and  that  for  the  linear  pressure  distribution  cases, 
pressure  vector  pf  will  have  half  the  number  of  elements  in  the  vectors  in  equation  (9)  required  for  the 
consistent  distribution  case  and  jT  f  ^T. 

Substituting  equation  (10)  into  equations  (5)  and  (8),  we  get 


P  ( x  )  = 


1 

271 


I 

J=1 


r  * 

J 

Aj 


,r  ikR 

^  f„  <V>dAJ 


pS 

.1 


-ikpc-  /  aT 

A.i  J 


T  ikR 

T  •  e-p-  dAj 
R  J 


(II) 


d 

'9Tk 


p(x‘) 


1 

2tt 


.1  =  1 


LAj 


T 
c  . 
“I 


An^t. 

k 


ikR 

(e-r)dAj 


P£  + 

“J 


* 

f  T  a  ikR 
-ikpc  J  a  3T-  (— )dAJ 
Aj  k 


02) 


Appendix  A  provides  the  details  of  numerical  integration  to  be  performed  in  evaluating  the  above 
integrals.  Collecting  the  terms  in  proper  order,  we  can  write  equations  (11)  and  (12)  in  matrix 
notation: 


(P  }  =  f H | { P  }  +  (C I { V  }  ,  (13) 

n  n  n 

In  the  above  equations,  n  indicates  that  we  are  dealing  only  with  the  components  of  pressure  and 
velocity  normal  to  the  surface. 

Let  [El  be  an  identity  matrix  of  the  size  of  [H];  then  we  can  rewrite  equation  (13)  in  convenient 
form: 

[E](P  }  =  [H]  (P  }  +  [G] (V  } 
n  n  n 

[[E]  -  [H]]{Pn>  =  [G]{Vn> 

[H](P  }  -  [G]{V  },  (14) 

n  n 

where 

[H]  =  [EJ-IHJ. 
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Equations  of  Motion  of  the  Structure 

The  equations  of  motion  of  a  structure  idealized  into  a  finite  element  model  are  written  in  matrix 
notation  as 


[M]  (x)  +  [  K  ]  {x}  =  {F  }  +  {fJ, 

e  r 


(15) 


where 

[M]  is  a  mass  matrix 
[K]  is  a  stiffness  matrix 
{ x}  is  the  vector  of  generalized  displacement 
{Fe}  is  the  vector  of  generalized  mechanical  forces 

{ Ff}  is  the  vector  of  generalized  fluid  pressure  forces  normal  to  the  surface  . 

In  practice,  stiffness  and  mass  matrices  of  a  structure  are  assembled  by  (1)  writing  element  stiffness 
and  mass  matrices  for  each  finite  element  in  their  local  coordinate  system,  (2)  transforming  all 
element  matrices  individually  in  the  global  coordinate  system,  and  (3)  arranging  and  assembling  all 
the  matrices  in  final  form. 

A  generalized  set  of  displacements  {u}  of  a  material  point  of  a  jth  finite  element  is  described  by  a 
matrix  of  basis  functions: 

{ u }  =  [aHU^}  , 


where 

{u}  is  the  (r  x  1)  vector  of  the  generalized  displacements,  r46 
[a]  is  the  (r  x  k)  rectangular  matrix  of  basis  functions 

{U  ^  }  is  the  vector  of  (k  x  1)  elements,  where  k  =  the  number  of  generalized  coordinates  defining 
the  finite  element. 

Using  the  energy  principle,  we  can  derive  the  element  mass  matrix: 

[Ml  =  f  P[a]T[ a]dV.  , 

where  3 

Vj  is  the  volume  of  the  finite  element 
e  is  the  mass  density  of  the  finite  element , 

The  mass  matrix,  as  defined  above,  is  the  consistent  mass  matrix;  however,  a  lumped  mass  model 
can  also  be  used.  The  difference  between  the  two  mass  matrices  is  significant.  In  the  consistent  mass 
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matrix,  the  deformed  configuration  used  is  the  same  as  the  one  used  to  derive  the  stiffness  matrix.  The 
lumped  mass  approach  uses  a  discontinuous  deformed  configuration  in  which  several  parts  of  a 
particular  finite  element  tend  to  vibrate  as  pistons,  each  with  different  amplitudes.  The  advantage  of 
the  consistent  mass  matrix  is  that  it  requires  the  structure  to  be  divided  into  fewer  finite  elements  than 
does  the  lumped  mass  approach.  The  disadvantage  is  that  the  consistent  mass  matrix  is  nondiagonal, 
and  most  of  the  eigenvalue  routines  available  accept  only  the  diagonal  mass  matrix.  However,  since 
the  problem  we  are  concerned  with  uses  the  inversion  routine,  our  use  of  the  consistent  mass  matrix  is 
not  a  disadvantage.  From  the  previous  discussions,  it  can  be  seen  that  the  proposed  plan  to  use 
consistent  velocity  and  pressure  distribution  to  solve  the  Helmholtz  surface  integral  is  a  simple  ex¬ 
tension  of  the  idea  which  led  to  the  development  of  the  consistent  mass  matrix. 

The  generalized  forces  can  be  obtained  using  the  virtual  work  principle.  Equating  the  total  virtual 
work  done  by  the  mechanical  forces  and  the  generalized  forces,  we  get  for  the  jth  finite  element: 

{ }T{6V^}  =  /  p(x.)<5V  (x.)ds.  , 

fn  j  •'s  -j  n  ~j  j 

where  Sj  is  the  surface  of  the  finite  element. 

Substituting  the  expression  for  p  and  Vn  from  equation  (10),  we  get 

T 

f  p(x.)6V  (x.)  dsj  =  {[  {P^}  [C.][a.]Tds.  }{6V^}  > 

Js  -J  n  “J  jls  J  J  J  j  .1 

j  j 

which  results  in 

<>-[/,.  h1":/  deil(V  '  'V'V- 

The  matrix  equation  of  the  generalized  forces  for  the  entire  structure  can  be  written  by  properly 
arranging  the  above  element  generalized  forces: 

{ Ffn }  =  [A]{Pn} .  (16) 

This  expression  gives  the  normal  generalized  forces  in  terms  of  the  normal  pressure  vector. 


As  explained  later  in  the  Solution  of  Coupled  Equations  section,  we  can  write 

{Pn}  =  [T]  [B]  {P} 
l F/n }  =  [T]  [B]  { Ff } . 

Substituting  equation  (17)  into  equation  (16),  we  get 

{Ff}  =  [B]T  [T]T[A][T] [B]{ P }  . 

Substituting  equation  (18)  into  equation  (15),  we  get 

[M] {x}  +  [K] {x}  =  {F  }  +  [B]T[T]T[A] [ T 1 [ B 1 { P }  . 

e 


(17) 

(18) 


(19) 


8 


•m 


Tk  5821 


Interaction  of  the  Fluid  and  Structure 


In  discussing  the  interaction  of  the  fluid  and  structure,  we  will  first  consider  a  simple  scattering 
problem. 


Scattering 


Before  we  deal  with  the  solution  of  equations  (14)  and  (19).  let  us  consider  the  case  w  here  a  train  of 
harmonic  plane  waves  traveling  along  the  axis  of  symmetry  impinges  on  the  structure.  Superscripts  i 
and  s  refer  to  the  incident  and  scattered  waves.  p'.v'  are  the  incident  pressure  and  velocity  vectors. 
ps,vs  are  the  scattered  pressure  and  velocity  vectors. 


Let 


j  (kx-Rot )  .  i  . 

=  v  e  l  =  v  i 

o  — x  ^x 


s  s 
v  =  v  n 


i  ,  s 
P  =  P  +  P 

vn  =  (v.ln)n  +  vS . 


(20) 

(21) 

(22) 


p,vn  are  the  total  pressure  and  velocity  vectors  normal  to  the  surface  of  the  structure.  _ix,  n  are  the 
unit  vectors,  respectively,  along  the  axis  of  symmetry  and  normal  to  the  surface. 


Scattered  components  of  pressure  and  velocity,  ps  and  vs,  satisfy  the  Helmholtz  surface  integral, 
yielding 


[ H ] { P  S }  =  [G]{vS}. 

Incident  pressure  and  velocity  are  related  by  the  equation  of  motion 

pf  =  3^  * 


(23) 


(24) 


For  harmonic  velocity  distribution,  it  can  be  shown  that 

pi(x,t)  =  -  cpfvi(x,t). 


Equation  (22)  relates  only  the  normal  components  of  velocity. 
Let 


i 

v  -n 


in 


n  n 
v  =  v  n  , 


(25) 
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where  vn  is  the  component  of  velocity  of  the  structure  normal  to  the  surface  of  the  structure.  Now  we 
can  rewrite  equation  (22)  as 


n 


v 


in 

v 


+  v 


(26) 


Substituting  equation  (26)  into  equation  (23),  we  get 

{pS}  =  [H]"1  [G]  {{v}n-{vin}}  .  (27) 


And  substituting  equations  (25)  and  (27)  into  equation  (21),  we  get 

{p}  =  -  cpf{vL}  +  [hJ1[G] (v}n  -  [H]1[G]{vln}  .  (28) 

The  analysis  of  an  incoming  plane  wave  at  an  angle  to  the  structure  is  given  in  appendix  B.  Ap¬ 
pendix  C  discusses  further  the  scattering  of  the  plane  wave  from  rigid  and  flexible  structures. 

Solution  of  Coupled  Equations 

We  are  dealing  with  acoustic  pressure  which,  by  definition,  is  a  small  variation  in  some  uniform 
pressure  field.  This  precludes  any  possibility  of  cavitation  in  the  fluid  near  the  structure,  and  we  can 
assume  that  the  velocity  of  the  structural  surface  is  the  same  as  the  velocity  of  the  fluid  at  that  point, 
and  that  the  pressure  exerted  on  the  structural  surface  is  the  same  as  the  pressure  in  the  fluid  at  that 
point.  It  is  this  condition  of  continuity  of  pressure  and  velocity  that  couples  equations  (14),  (19),  and 
(28). 


Before  we  couple  these  equations,  the  pressures  and  velocity  vectors  of  equations  (14)  and  (28)— 
which  consider  only  the  normal  component  of  velocities  of  the  grid  points  located  on  the  surface- 
will  have  to  be  transformed  to  the  generalized  force  and  velocity  vectors  of  the  entire  structure  in  the 
basic  coordinate  system.  The  set  of  grid  points  on  the  surface  of  the  structure  is  the  subset  of  all  the 
grid  poiints  on  the  structure.  Also,  the  pressures  and  velocities  in  equations  (14)  and  (23)  that  are 
normal  to  the  surface  of  the  structure  will  have  to  be  transformed  to  the  basic  coordinate  system. 

Let  [T]  be  the  transformation  matrix  relating  the  normal  variables  to  the  variables  in  the  basic 
coordinates.  Since  it  is  obtained  by  rotating  each  of  the  local  variables,  the  resulting  matrix  [T]  is 
orthogonal.  After  the  variables  of  this  subset  have  been  transformed  to  the  basic  coordinate  system, 
the  subset  must  be  related  to  the  total  set.  Conceptually,  this  is  performed  by  a  rectangular  matrix  [B]; 
for  example,  to  relate  the  vector  fx  1  to  vector  fx  ■>. 

R 

we  can  devise  the  matrix  shown  below.  J 1 

°  ■ 

Although  we  will  retain  the  matrix  [B]  (defined  above)  in  this  work,  in  practice  the  matrix  that  [B]  is 
to  be  expanded  to  (by  inserting  proper  rows  and  columns  of  zeros)  is  used. 

From  the  above,  then,  we  can  write 

(vn)  =  [T]{vS>  -  [T] [B] {v}  .  (29) 


10 


TR  5821 


Similarly, 

{P  }  =  [T] [B] {p}  .  (30) 

n 

Substituting  equations  (29)  and  (30)  into  equation  (14),  we  get 

[T] [ B ] { P }  =  [Hi1 [G] [T] [B] {v> . 

Since  [T]  is  orthogonal,  i.e.,  (T]1  =  [T]T,  the  above  equation  reduces  to 

(B]T[B]{p)  =  [B]T[T]T[hT1[G][T][B){v}  . 

[B]T[B]  is  a  diagonal  matrix  with  unity  along  the  diagonal,  except  for  the  points  (corresponding  to  the 
interior  points  of  the  structure)  on  which  there  is  no  fluid  pressure  and  which  are  zeros.  In  fact,  if  the 
above  equations  are  written  out  explicitly,  both  sides  will  be  identically  zero  for  the  interior  points.  In 
what  follows  we  will  replace  [B]T[B]{ P }  by  {P},  with  the  implicit  understanding  that  the  new  vector 
includes  the  pressures  at  the  surface  points  and  the  interior  points,  and  that  the  pressure  at  the  in¬ 
terior  points  is  zero. 

We  now  write 

{P}  =  [B]T[T]T[hT1[G][T][B]{v}.  (31) 

In  a  similar  way  we  can  transform  equation  (28): 

{P}  =  -cpf[B]TiB]{vi}  +  [B]T[T]T[HT1[G][T][B]{v}  - 

[B]T[T]T[Hl1[G][T][B]{vi}.  (32) 

It  is  understood  here  that  vector  { v'}  is  defined  in  the  basic  coordinate  system  of  the  structure  for  only 
the  surface  grid  points  and  that  {y‘}  is  the  vector’s  counterpart  for  the  total  set  of  grid  points.  Since  it 
can  easily  be  verified  that  equation  (31)  is  the  subset  of  equation  (32),  we  will  consider  the  general 
equation  (32). 

Substituting  equation  (32)  into  equation  (19),  we  get 

[M] (x)  +  [K]{x}  =  (Fe)  -  cpf[T1]T[A][T1]{v1}  + 

[T1]T[A][hT1[G][T1Hx}  -  [T1]T[A][Hl1[G]tT1Hvi}>  (33) 

where 

ITj]  -  [T]  [ B ) 

(x)  *  {v}  . 
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Direct  Analysis.  Since  we  are  considering  only  harmonic  forces, 

{F  }  =  {F*}eiWt 

e  e  (34) 

{x>  =  {X}eiWt 

{v.}  =  {v*}e'i'a,t. 

~ 1  “1 

Substituting  equation  (34)  into  equation  (33),  we  get 

[t  — 1  1  (35) 

-00?  [M]  +  [K]  -  iu»[T1]1[A][H]1[Gl[T1lJ{X}  = 

{ F*}  -  cpf  [T1]T[A]  [T11  {v*}  -  [T1]T[A][hI1[G][T1Hv*}  . 


The  solution  of  equation  (35)  yields  the  displacement  vector  {X},  which  can  be  substituted  back  into 
equation  (3 1 )  or  (32)  to  get  the  pressure  on  the  surface. 

Modal  Analysis.  The  matrix  to  be  inverted  in  equation  (35)  is  complex,  full,  and  sometimes  ill- 
conditioned  because  the  terms  of  the  stiffness  matrix  are  large  for  certain  ranges  of  driving 
frequencies.  Also  the  imaginary  part  of  the  matrix  is  small  compared  with  the  real  part.  An  alternative 
in  this  situation  is  to  employ  modal  transformation. 


Let  [  4>  ]  be  the  matrix  of  the  desired  N  modes,  and  let  these  modes  be  normalized  so  that 

[4>]T[M][<{>]  =  [E]  , 

where  [E]  is  the  identity  matrix.  We  know  then  that 


[<M  [K][<M  = 


X2 

1 


2  •. 


[MD. 


(36) 


(37) 


Let  us  use  the  transformation 


X2  J 

N 


(xl  =  [<J>]  {q>  ,  (38) 

where  {q}  is  the  vector  of  the  generalized  coordinates. 

Substituting  transformation  (38)  into  equation  (33)  and  premultiplying  the  resulting  equations  by 
[$]T,  we  get,  using  the  relations  (36)  and  (37), 

{q}  +  m°{q}  =  MT{Fe}  -  cpf[0]T[T1]T[A][T1]{vi}  + 

[<J>]T[T1]T[A3[HT1[GHT1][<|)]{q}  - 
[^>]T[T1]T[A][HT1[G][T1]{vi}  . 
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Since  the  forces  and  resulting  displacements  are  all  harmonic,  the  above  differential  equations  reduce 
to 


[[X]D  -gj2 [E]  -  ia)[4>]T[T1]T[A][Hl1[G][T1]  [<()]]  {q*}  = 
[$]T{F*}  -  cPf  [d)]T[T1]T(A][T1l  {v*}  - 
[<()]t[t1]t[a][hT1[g][t1]{v*}  , 


(39) 


where 

(q>  =  {q*}eiU)t 

(F  }  «  {F»]ela)t 
e  e 

(v  }  *  {v*}eia}t. 

~i  1 

The  solution  of  equation  (39)  will  be  simpler  than  that  of  equation  (35)  because  the  matrix  is 
smaller.  Using  equation  (38)  to  transform  the  solution  {q*}  back  to  {X},  we  get  the  displacements  of 
the  structure,  and  from  the  displacements  we  get  the  pressures  on  the  surface  of  the  structure. 

Examples.  Of  the  many  examples  solved,  the  solutions  of  two  problems  are  presented  here.  In  the 
first  problem,  as  shown  in  figure  3,  the  pressure  distribution  on  the  surface  of  the  rigid  sphere  is 
calculated  for  a  prescribed  velocity  distribution.  In  the  second  problem,  a  steel  sphere  submerged  in 
fluid  and  excited  by  a  point  load  is  considered. 

The  series  solution7  of  the  acoustic  pressure  on  the  surface  of  a  point-driven  sphere  has  been 
programmed.  Fifty  terms  of  the  series  were  retained  for  the  computation.  The  results  are  plotted  as  a 
series  solution  in  figure  4.  Two  different  finite  element  models  of  the  sphere  were  considered.  The 
model  used  in  conjunction  with  equation  (35)  represents  the  sphere  by  8  conical  shell  elements, 
whereas  the  model  used  with  equation  (34)  represents  the  sphere  by  18  conical  shell  elements. 

The  first  three  modes  of  the  18  element  model  were  retained  for  the  calculations  T  he  intention  was 
to  see  how  well  the  series  solution  compares  with  the  Fluid  Interacting  with  Sri  net u res  (FIST) 
program,  in  which  relatively  crude  finite  element  models  are  used.  The  agreement  among  the  three 
sets  is  excellent  everywhere  except  in  the  vicinity  of  the  point  load.  This  is  to  be  expected  since  even  for 
a  dynamic  response  of  the  sphere  in  the  vicinity  of  the  load  in  vacuo,  a  more  refined  finite  element 
model  than  the  one  used  here  is  required. 


Application  to  Computer  Programs 

The  computer  program  was  structured  to  solve  four  cases:  (1)  radiation  from  a  vibrating  cavity,  (2) 
scattering  of  a  plane  wave  traveling  along  the  axis  of  symmetry  from  a  rigid  cavity,  (3)  radiation  from 
a  vibrating  structure,  and  (4)  scattering  of  a  plane  wave  traveling  along  the  axis  of  symmetry  from  a 
flexible  structure. 
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Figure  3.  Pressure  Distribution  on  the  Surface  of  the  Sphere 
for  a  Given  Velocity  Distribution  ka  =  2 
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Also,  two  types  of  pressure  distribution  over  the  shell  along  the  generator  are  considered:  “con¬ 
sistent”  and  “linear.”  In  the  consistent  case,  the  pressure  is  described  by  a  third  order  polynomial.  In 
the  linear  case,  the  pressure  variation  is  linear  and  is  described  by  a  first  order  polynomial.  In  both 
pressure  distribution  cases,  the  velocity  variation  is  given  by  a  third  order  polynomial. 

The  geometric  shape  and  driving  frequency,  as  well  as  the  variation  of  velocity  and  pressure, 
determine  the  coefficients  of  the  (H)  and  [G]  matrices.  The  size  of  vector  { Pn}  in  the  case  of  consistent 
pressure  distribution  is  twice  that  in  the  linear  case.  In  the  consistent  case,  the  variables  are  the 
pressure  and  the  spatial  derivative  of  pressure  at  each  grid  circle.  In  case  (1),  for  a  given  velocity 
distribution  v„,  the  analysis  gives  the  pressure  distribution  for  both  the  consistent  and  linear  cases. 
For  the  cases  considered,  the  actual  magnitude  of  pressure  in  the  consistent  and  linear  cases  varies 
little,  but  since  only  a  few  segments  are  required  for  the  solution,  the  difference  is  larger  within  each 
segment,  as  shown  in  figure  3.  Also,  the  consistent  case  agrees  better  with  the  known  exact  solution.  It 
is  estimated  that,  for  satisfactory  results  in  case  (1),  the  length  of  each  segment  should  be  less  than  or 
equal  to  one-fourth  the  wavelength  of  the  velocity  distribution  along  the  generator  of  the  shell.  This 
observation  leads  us  to  believe  that  the  mesh  selected  for  finite  element  analysis  of  the  structure  will 
prove  adequate  for  the  solution  of  equation  (14). 

In  case  (2),  which  is  simply  the  extension  of  case  (1)  since  the  inclusion  is  rigid,  the  net  velocity  { v11} 
along  the  surface  is  zero.  This  shows  that  the  scattered  velocity  {vs}  is  simply  the  negative  of  the 
incident  velocity  {vin}.  The  fundamental  relation  (14)  gives  the  pressure  caused  by  a  scattered  wave. 
Therefore,  The  observations  for  case  (1)  apply  as  well  to  case  (2). 

Similarly,  case  (4)  is  an  extension  of  case  (3).  In  case  (3),  it  is  observed  that  for  the  solution  of 
equation  (35),  the  available  inversion  routine  does  not  give  a  reliable  inverted  matrix  for  problems  of 
a  size  larger  than  20.  In  this  case,  it  is  advisable  to  resort  to  the  solution  of  equation  (39).  The  results 
obtained  using  equation  (39)  for  a  spherical  shell  divided  into  eight  segments  and  driven  at  a  point 
compare  well  with  the  exact  solution  everywhere  over  the  sphere  except  in  the  vicinity  of  the  point 
load  (figure  4). 


Summary 

Analysis  of  scattering  and  radiation  of  acoustic  waves  from  an  axisymmetric  structure  in  an  infinite 
fluid  medium  has  been  presented.  It  has  been  shown  that  the  use  of  a  basis  function,  consistent  with 
the  one  used  for  axisymmetric  pressure  on  the  finite  surface,  yields  accurate  results  and  uses  fewer 
finite  surface  elements  than  does  piston-type  analysis. 
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Appendix  A 


Analysis  of  Axisymmetric  Shells 


Introduction 


Normally,  an  axisymmetric  shell  is  idealized  into  a  series  of  axisymmetric  conical  shell  finite 
elements.  The  finite  surface  elements  for  the  solution  of  the  Helmholtz  surface  integral  are  selected 
here  to  be  the  axisymmetric  conical  segments.  The  basis  function  that  defines  the  velocity  distribution 
on  a  conical  segment  is  selected  to  be  the  same  as  the  basis  function  used  for  axisymmetric  conical 
shell  finite  elements. 

In  this  section,  the  directional  derivatives  of  exp(ikR)/R,  which  are  used  in  equations  (11)  and  (12, 
and  their  numerical  solutions  are  dealt  with,  along  with  the  derivation  of  generalized  fluid  forces  and 
the  consistent  mass  matrix. 


Main  Analysis 


An  arbitrary  pressure  and  velocity  distribution  on  the  surface  of  the  cone  can  be  described  in  a 
series  form: 


p (s  ,<p)  =  )  p  (s)Cos  m<J>  +  /  p*(s)Sin  m<{) 
ta  ‘ — ,  m 

m=0  m=l 

00  00 

v(s,<{>)  =  /  v  (s)Cos  m<f>  +  /  v*(s)Sin  m<f> 
_  m  *- — -  ” 


(A-l) 


(A-2) 


m=0 


m=l 


We  will  consider  the  case  that  requires  the  pressure  distribution  corresponding  to  a  velocity 
distribution  given  by  only  one  term  of  the  infinite  series,  say,  vm(s)  cos  m$: 

p(s,<p)  =  p  (s)Cos  md>  (A-3) 

m 

v(s,<f>)  =  v  (s)Cos  m<j).  (A-4) 

m 

From  now  on  we  will  drop  the  subscript  m  of  pm(s)  and  vm(s),  and  we  will  evaluate  integral  equations 
(5)  and  (8)  for  an  axisymmetric  structure,  as  shown  in  figure  A-l: 


f  ikR  ikR 

P(x')  *  2ir J  ikpcv (x)  )dA  (A-5) 

'A 


TR  5821 


a  i  C  a2  1KK  a  1RK 

afj  (p(^1))  *  *fJa{p<2>  3^4  C^T)  •  ikpcv(i)  (iV)}dA 

R2  =  {  |x-x‘  |  •  Ix-x1 1 }  -  {(x-x1)2  +  (ySin<J>-y 1  Sin<J> 1 ) 2  + 

(yCoa<t>~y 1  Cos<P 1 ) 2  } 

=  {(x-x1)2  +  y2  +yl2  -  2yy' (SiniJiSinfjAcos^CostJ)1 ) 

=  {(x-x1)2  +  y2  +  yl2  -  2yy1Cos(<f>-<|>1)}.  (A-6) 

We  will  consider  that  point  x '  is  always  along  the  generator  of  the  shell,  which  implies  that  + '  =  0: 

R2  =  (x-x1)2  +  (y-y1)2  +  2yyJ (1-Cos +) .  (A-7) 

Define  the  normal  derivative: 


TR  5S:i 


Define  the  tangential  derivative: 


— —  =  v*t  =  (iv-+i-^rr+i  T-)  *  (_i  Cos  6  +  i  SinS) 
3t2  -  -2  Vjtr9r  — <J>r  3<t>  -x9x  -r  “X 


JAl 


=  -  Cos  B  +  Sin  B 

dr  dx 


Since  the  shape  is  axisymmetric  in  the  above  two  formulas,  we  can  substitute  y  for  r; 


=  Sin  9  +  Cos  0 

9n  9y  dx 


^ —  =  -  Cos  B  +  Sin  B  jr-  •  (A-9) 

9t2  9y  9x 

Substituting  equations  (A-3)  and  (A-4)  into  equations  (A-5)  and  (A-6),  dropping  the  subscript  m  of  p, 
v,  and  the  subscript  2  of  t2,  and  letting  $  '  =  0,  we  get 

.  L rV  o  ikR  ikR 

pfs1)  =  ^  J  J  (P(s)  Cosm<J>  ikpcv(s)  Cosm«t  }  yd<j>ds 


,2  ikR 


p, jj  (s 1 )  =  YU  J  J(p(s)  Cosm <p  3na  t~(~  R  )  "  ikpcv(s)  CosmiJ)  y^-(-e-~-)  yd<{)ds 


Assuming  that  the  total  surface  is  made  of  J  conical  finite  surface  elements,  we  get 


J  L 


p(s‘’  ■ w  Il(si*sj)p(sj)y(sj)dsj  +  /  I2(si’sj)v(sj)y{sj)dsjJ  (A- 10) 


J  L 


.,<*  -5J[f  I3(si,Sj)p(Sj)y(Sj)dSj  +/  IA(si»sj)v(sj)y(sJ)ds^  t(A-ll) 
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where 

Sj  refers  to  the  ith  grid  circle  where  pressure  is  to  be  found 

Sj  refers  to  the  field  point  on  the  generator  of  jth  conical  finite  surface  element  0<Sj<Lj 

y(Sj)  is  the  y  coordinate  of  the  field  point  Sj 

v(sj)  is  the  velocity  normal  to  the  surface  of  the  field  point  s^ . 

p>t  (s^  =  I  P^Sj)!  =  9t-jp(si)|  (Note  the  way  the  index  of  t  is  changed.) 


I 


i(s  , s  ) 


Cosmtj)  d<J> 


12(si’Sj) 


...  r  ikRj  * 

-  7_C  f  - )Cosm4>  d4> 

27T  —it  Rij 


I3<Si’Sj) 


aY 

2  IT  J 
-TT 


9n  3t . 
3  i 


ikRii 

(^— - )  Cosm<t>  dcf) 

Rij 


VVV  =  - 


ikpc 

2tt 


-7T 


j.  ikRii 

x—  (^— - -)  Cosing  d<j) 

3ti  Rij 


9 

3n. 

J 


Sin© 


j9y; 


+  Cos© .3 — 
j  9x . 
J 


9 


Cos0.3 —  + 
i9 


Sin0 


l9x . 


{(x^-x1)2  +  (yj-yi)2  +  2yiy^ (l-Cos^)}1^2 


A-4 


(A- 12) 


(A- 13) 

(A- 14) 


(A-15) 


(A- 16) 


(A- 17) 


(A- 18) 


TR  5821 


ji 


Polynomial  Evaluation 


Given  the  slopes  (if  and  the  ordinates  Yj  and  Yj+ ,,  and  if  the  curve  is  defined  by  a  third  order 
polynomial 


yj  =  C0+  Clyj  +  Vj  +  c3yj  > 


coefficients  c0  through  c3  can  be  determined  by  the  four  boundary  conditions  fij,  /if,  Yjt  YJ  + (See 
figure  A-2.)  The  resulting  expression  can  be  rearranged  so  that 


yj  »  (l-3x^+2xJ)Yj  +  Lj(x-2x^+xJ)$j  +  (3xZ-2xJ)Yj+1  +  (-x2+x3)6j+1 ,  (A- 19) 


where  3 

X  „ 


j 


Rewritten, 


y.  =  A1  (1)Y  +  A1  (2)Bj  +  A2j(l)Yj+1  +  A2j(2)Bj+1  ,  (A-20) 


where 


Alj (1)  =  l-3x2+2x3  A2j(1)  =  3x2"2x3 


Alj (2) =(x-2x2+x3)Lj  A2j(2)  =  (-x2+x3)L  • 


(A-21) 


where 


yj  -  v5jc°s'j +  vini  ■ 


yj  =  A1  (2)B‘  +  A2j(2)02 


(A-22) 


(A-23) 


A- 5 


,y ...  »  ..ST: 


TP  S  J  1 


6.  =  y.  =  (1-4x+3x2)6!  +  (-2x+3x2)B2 
J  dSj  J  J  J 


e  =  90  +  a:2  +  tan1^ .  . 

J  J  J 


Velocity  Distribution 


v(s.)  =  A1.(1)W.  +  A1 . (2)a .  +  A2 . (1)W  .  +  A2.(2)«  . 

J  J  J  J  J  J  J+l  J  J+1 

=  {ai.}t{dJ  +  (a2 . }t{d  }  , 

33  33+1 

where  W,  =  displacement  normal  to  the  surface  at  jth  grid  circle. 


C  =  -i-  W  {Al.}x  =  (Al.(l),  Al.(2)} 

3  dSj  3  3  3  3 


W.  =  -j-  W. 
3  dt  3 


{A2  }A  =  {A2.(l),  A2j  (2) } 


(d.) 

J 


W. 

.3 


Pressure  Distribution 

p(s.)  =  A1.(1)P.  +  Alj  <2)Pj  +  A2ja)P.+1  +  A2.(2)P^+1 

=  {A1.}T{S.}  +  {A2.}T(S.+1}, 

3  j  3  3+1  ’ 

where  Pj  =  pressure  normal  to  the  surface  at  jth  grid  circle. 

p1  =  _d_  p 

3  dSj  j 

If  only  the  linear  variation  is  to  be  used, 

p(s . )  =  cl .P  .  +  c2  .P  .  .  , 

3  33  33+1 


where 


s . 


s  . 


cl.  =  1 - T*4-  c2.  =  -r-1 

J  L.  J  L. 

J  3  3 


(A-24) 

(A-25) 

(A-26) 


(A-27) 

(A-28) 

(A-29) 

A- 7 
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Evaluation  of  I,(s,«  Sj) 

Evaluation  of  Directional  Derivatives 


;ikRij 

Rij 


CoskRij  +  SinkRij 
Rij 


»  ikRij  ,  ikRij  . . 

3  .e  J  ..  _  3  ,e  ^Rij 

3n . '  Rij  '  ’  O’  ’  ;  = 

J  J 


3  Rij  Rij  3n. 


3  •  elkRlJ 

—  - __) 

3t.  Rij 
1  J 


3  /eikRij,3Rij 

i  •  •  N  r*  5  •  )  * 


3Ri j  Rij  3t. 


3RTT<£^rj~)  =  RlJz(Rij-iReikRij  -  eikRij> 


From  equations  (A- 16)  and  (A- 17), 


e ikRij  eikRiJ 

*  ik - —r~. —  -  - —  :-r? 

RlJ  RlJ* 


_l_8ij  =  _l_Rij  .  Cos8j  +  Sin0j- 

j  j  j 


+  Sinei-^j- 


.iRil  +  Sine.^il 


£ 

- — Rij=  -  Cos 9  . +  Sin9.:» 
3t.  i  3y.  l  3x 

i  j 


Rij2  =  (Xj-x^)2  +  (yj-y^)2  +  2y^y^  ( l-Cos<(>) 


2Ri  j 


3x . 
J 


2(x . -x .  ) 
J  i 


9Ri  i 
3x . 
J 


x .  -x . 

_J _ L. 

Rij 


A-8 


(A-30) 
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x.-x. 

_i _ I 

Rij 


3Rii  -ir  „ 
ay.  ‘  ^ij{yj-yiCos0) 

=  Rrj{yryjCos^} 


3  /e*kR*\  x;_x^  y-  y.Cos4> 

30j(  Rij  }  "  (ifeRiJ-D-grp-fCoaejC-^-)  +  Sin9  j  (— J-  ~ - 


)}  (A-31) 


ikRi  j 


3  /ei,w'VJ  s  . eikRiJ  .  y.-y.Cosj 


3t.(  Rij  )  (ikRij_1)  ^Tij7  {'Cos0i  -Srfj—  +  Sin9i  .  (A-32) 

Simplifying  the  terms  shown  below,  we  get 


ikRi  j 

(ikRij-1)  _  -^TjjHCoskRij  +  iSinkRij 


ii) 


“(  RipCoskRij'  RXjSinkRij)  +  i(~CoskRij  -  SinkRij) 

=  -(  CpskRij+kRijSinkRij,  .  . ,kRi jCoskRi i -SinkRi i „ 

1  £iF^  >  +  l(  Rij*J - (A-33) 

snufi7  ■  BiF)e‘kRlJ,-<ik<sF  •  bf1  * (-  Mr  *  sfpe>*ikRii 


j^2  ^  3k 

Rij?  +  RiJT  _i(Rlp)>  ( Co  skRi  j+i  SinkRi  j  ) 
{(2Z|TJ^“)COSkRij+|rpSinkRij}  +  i{(--R'lj4j2)SinkRij  -|jpC0SkRij} 


^  ikRij  r  2  Z 
8nj3t^  ^  Rij  ^  ^  RijA^  )CoskRij  +  SinkRij} 


2„.  .2 


3-k2Ri12 

^  Rij *  SinkRij  -  ^|yj  CoskRij}|  {-(y^y^Cos^OCosSi  + 


A-9 


i 
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(x  ,-x  .  )Sin0i}  •  — tt  {(x.-x,)  Cos9j  +  (y.-y.  Cos  <p)  SinOj}.  (A-34) 

1  j  Rxj  j  1  31 

Substituting  equations  (A-30)  through  (A-34)  into  equations  (A-12)  through  (A-15*  and  separating 
the  real  and  imaginary  parts,  we  can  write 


where 


Il(si’sj)  =  Ill(si*Sj)+iI12(si*Sj) 

(A-35) 

I2(si’Sj)  =  I21(si’Sj)+iI22(Ri’Sj) 

(A-36) 

I3(si’Sj)  =  I31(si’Sj)+iI32(si’Sj) 

(A-37) 

I4(si’Sj)  =  I41(si’Sj)+iI42(si’Sj)' 

(A-38) 

CoskRlj^iSinkRij}|-(^ ^ Cos0 j+(y  _y  Cos()))  Sin0 jj 

Cosm<}>d<J>  (A-39) 

TT 

I  fs  si-1  f  \  kRUC?.s.kR,y-StnkRij  > 

i12l  i’  j  2tt  J  *  RijJ  • 

J  -TV 

•  {  (x  ,-x  .  )Cos0j+(y  ,-y  .Cos<j>)  SinOj  }Cosra<J>d<j> 

3  i  J  i 

(A-40) 

— TT 

(A-41) 

*22  ’  *  ^ 

-TT 

(A-42) 

131  -  j-  )'{ 2)CoskBij^RUSJt,k«U,  „  f2  (AJ3) 

— TT '  J  ' 

132  .  X  >SinkRlj-3kRlJCoskRlJ  }  CM^  (A-44, 


A- 10 


TR  ss:i 


FI 


(x^-x^ )  Sin9i-(yi-yjCos4))Cos9i 


F2 


(x  j-x  .)Cos9j  +  (yj-yiCosc^)Sin9j 


I 


41 


TT 

kpc  /»  (kRijCoskRij  -SinkRij 

2"  J  l  Rip 


{ 


(x^-Xj )  Sin0i-(yi-y^.Cos<{))Cos9l|Cosm4>d<() 


(A-45) 


T  _  kpc  r  t CoskRi j+kRij SinkRij ^ 
42  ~  2n  _7T<  Rij3  f 


■|  (xi-x^.)Sin0i-(yi-y_.Cos<J))Cos0i|cosin({)d^  . 


(A-46) 


Assumption  I.  Consistent  Pressure  Distribution 

Substituting  the  integrals  I,  through  14  and  equations  (A-26)  and  (A-27)  into  equations  (A- 10)  and 
(A-ll),  we  get 


P(s.)  =£[{G1(i.j)}T{D  H{G2(i,j) 
j-1  J 


T 

>  {Vi} 


+{Hl(i,j)}T{Sj }+{H2(i,j)}T{Sj+1}] 


(A-47.1) 


P,t(s.)  =  £  fc3(i,j)}T{D.  }  +{G4(i,j)  }T{Dj+^ 


+  {H3(i,j)}T{Sj)+{H4(i,j)}T{Sj+^  . 


(A-47.2) 


The  two  preceding  equations  (A-47)  can  be  combined  in  one  matrix  equation: 


A-ll 
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where  i  =  1,2 


J 

{s . >=  22 

1  j=i 

|"  {Gl(i,  j) }T  {G2(i,j)}T 

,<V  | 

^{G3(i,j)}T  {G4(i,j)}T 

{Vi’j 

+ 

(HI (i , j ) }T  (H2(i,j)}T 

{ H3 ( i , j ) )T  (H4(i,j)}T 

'V 

+  1) 


(A-48) 


i  =  1,  2,  .  .  .(J+l) 

{Gl(i,j)}T  =/  I  (s  ,s  )y(s  ){Al.}Tds, 

j  J  j 

(A-49.1) 

L. 

{G2(i, j)}T  ^  I2(s1,sj)y(sj){A2j}Tdsj 

(A-49.2) 

ri 

(G3(i, j) }T  =J  I4 (s^  >sj )y (s^ ) {Alj }^dSj 

(A-49.3) 

L. 

J 

(G4(i ,  j )  }T  =  f  I  (s  ,s  )y(s  ){A2.}Tds 

J  0  ^  1  J  J  J  J 

(A-49.4) 

L. 

{Hl(i,j)}T  Il<Si*Sj >y(Sj){Alj }Tds^ 

(A-49.5) 

Lj 

[H2(i, j) }T  I1(s1,sj)y(sj){A2j}Tdsj 

(A-49.6) 

(H3(i,  j)  }T  =J  I3(si.sj)y(sj){Alj}Tdsj 

(A-49.7) 

Lj 

(H4(i,j)}T  -/  I  (s  ,s  )y(s  ){A2  }Tds  . 

n  J  1  J  J  J  J 

(A-49.8) 

A-12 
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Equations  (A-48)  represent  2(J  +  1 )  equations,  which  are  written  in  matrix  form  as 

[E] {s}  -  [G]{D}+[H]{S},  (A-50) 


where 

[E]  is  the  unit  matrix,  (2 J  +  2)  size 

[G],  [H]  are  complex  square  matrices,  (2J  +  2)  size 

{S},  {D}  are  the  vectors  of  generalized  pressure  and  velocity,  respectively,  of  (2J  +  2)  size. 
Assumption  2.  Linear  Pressure  Distribution 

Substituting  the  integrals  I,  and  I2  and  equations  (A-26)  and  (A-28)  into  equations  (A- 10),  we  get 


j 

P(s.)  ^{Gl(i,j)}T{Dj}+{G2(i,j)}T{Dj+1}+Hl(i,j)Pj+H2(i,j)Pj+1J  ,  (A-51) 


where 


ii  i 

Hl(i,j)  -J  I1(si,sj)y(sj)cljdsj 


H2 


(i.j)  =J"  I1(si,sj)y(sj)c2jdsj  > 


(A-52.1) 


(A-52.2) 


where 


cl  =  i  _  -J-  c2  =  -j 
j  1  L4  ’  CZj  L. 

J  j  J  3 


The  total  of  (J  +  1)  equations  (A-51)  are  written  in  matrix  form  as 

[E] (P)  =  [G]{D}+[H](P},  (A-53) 

where 

[E]  is  a  unit  square  matrix  of  (J  +  1  ]  size 

[G]  is  the  complex  matrix  of  (J  +  1),  (2J  +  2)  size 

[H]  is  the  square  complex  matrix  of  ( J  +  1 )  size 
{P}  is  the  vector  of  pressures,  (J  +  l)size 

{D}  is  the  vector  of  generalized  velocity  of  (2J  +  2)  size  . 


A- 13 
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Numerical  Integration 


It  is  apparent  that  the  integrals  (A-40)  through  (A-47),  (A-49*i(i  =  1,8)),  and  (A-52»i(i  =  1,2))  can 
be  evaluated  only  by  using  some  numerical  scheme.  Here  the  Gaussian  integration  formula  is  used 
because  the  integrands  are  analytical  expressions,  which  makes  it  easy  to  calculate  their  values  at  any 
required  point,  and  beccause  the  number  of  calculations  to  be  performed  for  a  required  accuracy  are 
fewer.  The  Gaussian  integration  method  reduces  the  integral  to  the  summation  of  the  weighted  values 
of  the  integral  at  the  preassigned  number  of  points.  The  Gauss-Legendre  integration  formula  is  ex¬ 
pressed  as 

In 


f  y(x)dx 

i= 


Aiy(xi). 


(A-54) 


Points  Xj  and  the  corresponding  weighting  function  Aj  are  given  in  standard  references  on  numerical 
analysis.  The  integrals  (A-49)  and  (A-52)  are  not  in  the  form  of  equation  (A-54);  as  a  result,  the  in¬ 
tegrals  require  transformation.  The  above  integrals  are  in  the  form 


Let 


L 


s  =  ^(1+s1), 


(A-55) 


which  transforms  the  above  integral  to 


L  1  n 

f  f(s)ds  =  J  f  f(t(l+s1))ds1  f(W  )). 

“0  Z  i=l  1 


(A-56) 


Since  f(Rjj($))  and  cos  m$  are  both  symmetric,  integrals  (A-40)  through  (A-47)  have  the  following 
property 


V 

j~  f  (Rij  (4>)  )Cosm4>d<t>  =  2j  f  (Rij  (<£)  )Cosm<j>d<l> 


In  order  to  use  the  Gauss-Legendre  integration  formula  for  an  alternating  integral  like  the  above, 
we  should  first  reduce  the  integral  to  the  sum  of  the  integrals  over  each  wavelength  and  then  use  the  n- 
point  Gauss-Legendre  integration  formula  over  each  of  the  subregions: 


Let 


(Rij  )Cosm<}>d<j>  = 


n 


’E 

4=1 


7T  i 
m 


f  (Rij  )Cosm<j>d<J>  . 


“(i-1) 

m 


♦  i+I/2)2,  <J>  -  ^0jH-2i-l)  . 


A-14 


A- 15 
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Generalized  Forces 

The  acoustic  pressure  distributed  over  each  finite  surface  element  has  to  be  expressed  as  the  joint 
forces  {Ffn},  as  shown  in  equation  (16). 

Total  virtual  work  equals  /n/L  D  6u  yd*ds 

«HfT  O  K  *  # 


=“  (Q1"1  5v/?  +  Q2m6  *oc™  +Q 1 m  sif*  +  Q2">6;  m}  + 

.  °  1  J  1  J  2  j  +1  2  J  +1 


?{Qim*  SV^+Q2™*  6 "  crm*+Qlm*  6y"*  +  Q2«*a‘»* 
O  1  :  1  .  2  ■  2 

J  J  J+l  2  J+l 


}  • 


(A-58) 


Substituting  equations  (A-l)  and  (A-2)  into  equation  (A-58),  using  the  relations  of  equation  (A-26) 
and 

f1’  Cosmiji  Cosn<|>d<ti  =  errfimn 

-it 

Sinmit>  Sinn<J>d<J)  =  ETrdmn 

-TV 


Sinm<{>  Cosn4>  d<(!  =  C 

-TV 

,  r5  in  5  0 

where  e  =  ,  n 

1  m  i  0 

and  then  comparing  the  coefficients  of  the  like  virtual  generalized  displacements,  we  get  the  ex¬ 
pressions  for  the  generalized  forces,  as  shown  in  figure  A-3:  (Note  that  superscript  m  is  dropped  and  j 
is  added.) 


Q^=  evr/p  (s)y(s)Al  .(l)ds 
1  om  J 

q2Js  eir/p  (s)y(s)Al .  (2)ds 

1  „  m  J 


Q1 j*  Ew/p  (s)y(s)A2 .(l)ds 
2  o  m  J 

2  4  T 

Q  J=  Eir/p  ( s)y(  s)A2  .  (2)ds  .  (A-59) 

2  o  ^ 


A-16 


.zESH 
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In  the  above  equations,  the  only  assumptions  we  have  made  concern  the  velocity  distribution,  given 
by  equation  (A-26),  and  the  definition  of  the  positive  direction  of  the  velocity  and  pressure  (along  the 
outer  normal  to  the  surface).  Although  the  generator  of  the  axisymmetric  segment  is  shown  to  be 
straight,  equations  (A-59)  admit  any  variation  of  y,  straight  or  curved.  No  restriction  on  the  pressure 
variation  of  the  axisymmetric  finite  surface  is  imposed  in  equations  (A-59). 

Depending  upon  whether  the  pressure  distribution  is  consistent  (equation  (A-27))  or  linear 
(equation  (A-28)),  one  can  analyze  the  integrals  of  equations  (A-59)  exactly  if  the  analytic  expression 
for  y(s)  is  available  (e.g.,  the  generator  is  straight)  or  numerically  if  y(s)  is  arbitrary.  In  the  analyses 
that  follows,  the  generator  of  the  segment  is  assumed  to  be  straight,  which  means 

y(s)  =  C1.-Y.+C2.-Y.  ,  .  (A-60) 

J  J  J  J+1 


Cl ,,  C2j  are  given  by  equation  (A-29). 


Assumption  I.  Consistent  Pressure  Distribution 
As  defined  previously, 

P  (s)  =  A1 . ( 1 )P ,+Al . ( 2)p!+A2 . ( 1 )P .  +A2  (2)P1 

m  J  J  J  J  J  ,1+1  j+1 

=  f  A1 .  >T{S  . }  +  { A2  .  }T{  S  .  (A-27) 

J  J  J  J+1 


Substituting  equations  (A-60)  and  (A-27)  into  the  integrals  of  equation  (A-59),  we  get 


(A-61) 


As  explained  in  equation  (17),  the  above  matrix  for  the  jth  element  can  be  transformed  into  the  basic 
coordinate  system,  so  that  the  total  matrix  can  be  built  from  these  element  matrices.  Here,  however, 
an  alternative  approach  is  taken  to  write  explicitly  the  total  matrix.  At  the  jth  node  point,  we  can  write 
the  following  relation: 


A- 18 


TR  ss:i 


Q  .  =  -Iq^Sin'oc  .+  Q1  (  j  1  ^SirT*  .  I 
XJ  \  1  J  2  J  -  If 

Qyj'K’i  COS’“  i  *  Q2<i"1>C“‘“i-,} 


■V  ^  *  «2<J"n 


Substituting  equation  (A-61)  into  (A-62)  and  rearranging  the  terms,  we  can  write 

Qy j  =  A( 3 j  —  2 ,  2j-3)Pj_1+A(3j-2,2j-2)p|_1+A(3 j-2,2j-l)P. 
+A(3j-2,2j)Pl+A(3j-2,2i+l)P.+,A(3j-2 ,2i+2)p!  , 

j  j  1  J  j  + 1 

Q  .  =  A( 3 j- 1 , 2 j -3 )P  +A( 3 j  —  1 , 2 j -2 )p!  ,+A( 3 j -1 , 2 j -1 )P . 
x3  j-1  1_1  '  '  J 

+A( 3  j  —  1 , 2 j )pj  +A(3j-l,2i+l)P.+  1+A(3j-l ,2j+2)pj  +  1 


Q3.  =  A(3j,2j-3)P._1+A(3j,2j-2)Pj!_1+A(3j  ,2i-l)P 


+A(  3  j  ,  2  j  )P  j  +A(3j,2j  +  l)P.  +  1+A(3.i,2j+2)r!  +  1 


9_ 

140  “j- 


A(  3 j~2 , 2 j -3 )  =  ~  ^ +  Y.) 


A(3j-2,2j-2)  -  L-l  *  ^ 


\  70  60 


A(  3  j  —  2 , 2  j  —  1 )  =  —■  Z  Y.  +  fv./Z.  +  Z.\  +  Jr-  Z.Y. 

35  j-1  j-1  7  j-1  j)  35  j  j+1 

A(3j-2,2j)  -U  jZ  Al  ♦  L.V  LjZ  /!i  !j±l 
J  J  \  60  28/  3  3  \  28  60 

A(3j-2,2j+l)  "  I45  Z.(Y.  Yj+iJ 

A(  3j  —  2  ,  2j+2)=  ~t,.Z./Yj  +  Yi-t-l) 

J  J'60  70  / 


(A-62) 


(A-63) 
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A(3j-1,  2 j— 3)  = 
A(3j^l,  2 j-2 )  = 
A( 3 j  —  1 ,  2 i-1)  = 

A( 3 j — 1 ,  2 j  ) 

A( 3 j-1 ,  2 j+1 )  = 
A( 3j-l ,  2 j  +  2 )  = 

A( 3 j , 2 j — 3 ) 

A(  3  j  ,  2  j — 2  ) 

A( 3 j , 2 j-1 ) 

A(  3 j , 2 j ) 
A(3j,2j+1)  - 

A(3j,2j+2) 

where  R.  =  erL, 
J  : 

Z.  =  e*L' 
J  : 


9 

140 


iY- 

35  j-1  j- 


-L.  R. 

j-1  j-1 


\  3 

'  ‘  35  RjYj+l 


Vi(Yi  *  Yi) 

Yirj  J'hzi ♦  h) 

3  J  y  70  60/ 

*  i  ■  -,)  • 

( '■  *  ’-) 

+  !i±i\ 

J  J\6°  70  ) 

2  /  Y .  Y.  \ 

i-1  f  -LZl  +  _I  ]  C71 
J  l  60  70  J 

?  1  /y.  ,  +  Y  A  £TT 

SHJ'1  v 

(h=l  +  +  Lj  (li  +  li±l)  } 

J  \60  28/  3  \28  60  / 


2 _ D 

140  J 


-L3 
280 


l?  Ylti  *  Ll\  -  *  liiiV 

3  \280  168 J  J\168  280  J 

](h  ♦  !m) 

J\70  60  / 


C7T 


L3 


280  \ 
2 

Sin". 

I  J 

.  Cos “? 


GTT 


,Yj +  Vi  r 


(A-64) 
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Writing  all  the  generalized  forces  Qxj,  Qyj,  Q3j  (j  =  1-J)  as  a  vector,  and  using  the  relations  of 
equation  (A-63),  we  can  define  the  following  relation: 

(Ff)  -  [*]  {s},  (A-65) 

where 


Matrix  [A]  3J  x  2J  is  highly  banded  and  its  coefficients  are  defined  by  equation  (A-64). 

Assumption  2.  Linear  Pressure  Distribution 
As  defined  previously, 

p  (s)  =  Cl.P.  +  C2.P.  ,  .  (A-28) 

ra  J  J  j  j+1 

Substituting  equations  (A-60)  and  (A-28)  into  the  integrals  of  equation  (A-59),  we  get 
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Substituting  equation  (A-66)  into  equation  (A-62)  and  rearranging  the  terms. 


we  can  write 


Q yj  =  A(3j-2,j-l)  p  1+A(3j-2,j)P.+A(3j-2,j+l)P. 

J  1  J  j+1 

Qxj  =  A( 3j-l , j-1 )P  +A(3j-l>j)p.4A(3j-1,j+I)p. 

J  1  J  j+1 

Q3j  =  A(3j,j-l)P._1+A(3j,j)Pj  +A(3j,j+l)P.+l 


(A-67) 


/  Y  V . 

A(  3 j-2 , j-1 )  =  Z  (  _±±  +  _1 

J  1  V  15  12 


A(3j-2,  j)  =  Z  H_±±  *  U  \  +  Z./4_Y.  +  liil 
J  1  V  12  15  V  J  15  J  12 


A(3j-2,j+l)=  Z.  +  ]j  +  I  ) 

J  \i2  is  y 

4<y-I’J-1>-  v>  l#1  *  y 

■  Vi  Yi)  *  Ri 

A(3j-1,  j+1)*  R.  +  ]j  +  l  ] 

J  l  12  15  y 

AOi.M)  — 

A(3j  ,  j)  =  -  /^-l  +  ewL*  +  Ii+l\ 

J  1  l  60  30  '  J  \  30  60  / 


30  60 


A(3j,j+1)  = 


e"Lj  (!i  -  !i±i\ 


60  60 


.  2 

-eirL.  Sin  n . 
J  J 


eirL .  Cos  <*  .  . 
J  J 


(A-68) 
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In  the  matrix-vector  notations,  the  generalized  forces  given  by  equation  (A-67)  are  written  as 

{ Ff }  =  [aJ  {P}  ,  (A-69) 

where 


r  i 

r 

Qxl 

pi 

Qyi 

P2 

Q31 

II 

Q3j 

P. 

J 

L  < 

Matrix  [A]  3J  x  i  is  highly  banded  and  its  coefficients  are  given  by  equation  (A-68). 


A-23/A-24 
Reverse  Blank 
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Appendix  B 

Analysis  of  an  Incoming  Plane  Wave 
at  an  Angle  to  the  Structure 

When  a  plane  wave  impinges  on  a  surface,  the  component  of  the  incident  plane  wave  norma'  to  the 
surface  scatters  off  the  surface.  Simultaneously,  because  of  this  acoustic  pressure,  the  structu.",  if 
flexible,  vibrates  and  radiates  acoustic  waves.  (If  the  structure  is  rigid,  this  vibrating  radiation 
component  is  absent.)  The  net  acoustic  pressure  is  the  sum  of  the  incident,  scattered,  and  radiated 
components.  Often  the  scattered  and  radiated  components  are  lumped  together  in  the  total  scattered 
wave. 

Distribution  of  the  incident  velocity  of  a  plane  wave  normal  to  the  surface  along  the  circumference 
of  a  circle  is  expressed  in  a  Fourier  series.  We  obtain  the  total  scattered  pressure  due  to  this  incident 
velocity  distribution  and  then  add  all  the  components  (incident,  scattered,  and  radiated)  to  get  the 
total  acoustic  pressure. 

The  particle  velocity  of  a  plane  wave  traveling  along  the  negative  y,  direction  is  given  by 

ll  -  iy  Vo  eK^rt)  ,  (B-l) 

where iy)  is  the  unit  vector  along  the  y,  -  axis.  We  will  drop  e~ja,t  from  the  analysis  that  follows. 

We  can  write  the  transformation  of  y,  into  x-y  coordinates  as 

y  *  (y-yQ)  Cos*  -  ( x-xq)  Sin*.  (B-2) 

As  shown  in  figure  B-l , 

y  ■»  r  Cos$.  (B-3) 


Figure  B-l .  Definition  of  Incident  Plane  Wave  With  Respect  to  Axisymmetric  Structures 


B-l 


mm 
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Substituting  equation  (B-3)  into  equation  (B-2),  we  can  write 


jky  _  jk(x  Sin*-y  Cos*)  #  jk(rCos«Cos<)>) 


e  1 


-jk  x  Sin*. 

•e 


(B-4) 


For  a  given  x0,y0,  <*, let 


c  =  ejk(xoSin«-yoCos«)  . 


(B-5) 


Also,  we  can  expand1  the  second  terms  of  equation  (B-4)  into  a  Fourier  Bessel  series: 


ellr?os.  Cos*  _  2  j  ^Cos*)  Cosmt _ 

m=o 


(B-6) 


where 


1  if  m  =  o 
e  =  . 

m  2  if  b  f  o  . 


Substituting  equations  (B-5)  and  (B-6)  into  equation  (B-4)  and  then  into  equation  (B-l),  we  get 


v1  =  iyVo.C.  je-jkxSin“-f  emjmJm(krCos=)  Cosing. 
~  (  m=o 


(B-7) 


We  will  derive  the  expression  for  Vin,  which  is  the  component  of  incident  velocity  normal  to  the 
surface: 

ryin  ..I  .  , 

V  =V  •n=Ai.n 


where  A  is  given  by  equation  (B-7)  to  be 

-jkxSin*  .n 


A  =  V  C)e 
o 


I  e  j  J  (krCos*)  CosmA 
m*o  mJ  m 


We  can  express  the  unit  vector  iy  and  the  unit  normal  to  the  surface  n  in  x-y  basis  as 


(B-8) 


which  gives 


i  =  Cos*i  -  Sin«i  „ 
-y  -y  -x 


n  *  -SinQi^  +  Cos0Cos$i  +  CosOSin^i^  , 


i  • n  ■  Sin*  Sin©  +  Cos0  Cos«Cos$ 


(B-9) 

(B-10) 

(B-l  1) 


Substituting  equation  (B-l  1)  into  the  above  expression  of  Vin,  we  get 


,in 


A{  Sin* Sin©  +  Cos*Co80Cos<J> } , 


(B-l  2) 


where  A  is  given  by  equation  (B-8). 


B-2 
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In  the  following,  we  will  try  to  simplify  the  expression  (B-l  2): 


( 

V  \ 


e  ^kxSin  Z  e  jmJ  (krCos*)  Cosm tpl 
m=o  mm  j 


{Sin*  Sin©  +  Cos*  Cos©  Cos<J>} 

=  Sin®  Sin©  e  ^kxSin  m|0  e  mjmJm(krCos*)  Cosm<)) 

+  Cos*  Cos©  e  ^kXSlnmf0  em  jmJm(krCos«)-|{Cos(m+l)4>  +  Cos(m-l)  $} . 
Rearranging  the  terms,  we  can  rewrite  the  above  as 

vljl  =  l=o  +  Bm(r,x,0)Cosm<(.  ->•  Am(r,x,9)Cos(m+lH} 


V  C 

Am(r,x,0)  =  -j- 


CosQCos*  e  J  kxS in* ^rCos* ) 

B  (r,x,0)  =  V  C  Sin© Sin*  e_JkxSin<rj  (krCos*)  , 

m  o  m 


We  rewrite  equation  (B-l  3)  as 
..in 


but 


■  Cos©  and  ~  *  Sin© 

~  -  Cos©  r^-  +  Sin©r^  • 
ds  3x  3r 


(B-l  3) 


(B-l  4) 


(B- 15) 


(Bo  +  liAl)  +Jl  Jj*  jVl  '  Am+1  +  jBm|  Cosm +  * 

Next  we  derive  the  expression  for  d/(ds)Vin  by  differentiating  equation  (B-l 3): 

d^ln  *m“o  emjm  j  ds(A®^  <Cos(m-l)*  +  Cos(m+l)*)  +  £(BTO)  Cosing ,  (B-16) 

By  the  chain  rule, 

d  ^  3x  ^  3  3r  t  3 

ds  3s  *  9x  3s  3r  ’ 


We  will  use  the  following  relations: 


|e  jl«Sin«j  >krCos*)(  -  kCos*  e_JkxSin<rj  (i^cos*) 


3xp  'j 

3“  | e“ J kxSln“ (krCo8«)|  -  kCos*  e"^kxSin*J^(krCos«) , 


B-3 
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where 


Ji  =  7  j  Jm-l(krCo8tt)  "  Jnfl(krCo8<r)S 


(B-17) 


1  d  .  o  ,  _  . .  -jkxSin* 

A  =  T-  A  =  — —  k  Cos0Cos«  e  J 
m  ds  m  2 


)-1Sin*Cos0J  +  Cos“Sin0J^( 

j  J  m 


B1  *  -1  b  *  V  Ck  SlnGS in*  e"jkx  Sln<1 
m  as  m  o 


-iSin*Cos0J  +  Cos*Sin0J 
J  m  u 


(B-18a) 


(B-18b) 


We  rewrite  equation  (B-16)  as 


vm 

ds 


K +  2i*lJ  2)"'x  jAi-i  -  >Li  *  rtf 


Case  I.  Wave  Traveling  Along  Axis  of  Symmetry 


(B-19) 


The  expression  of  the  velocity  Vin  and  the  rate  of  change  of  velocity  V'in  can  be  obtained  from 
equations  (B-15)  and  (B-19)  by  substituting  =  -90°.  which  will  align  the  y,  -  axis  with  the  axis  of 
symmetry: 


A  =  o 
m 


B«a  O  If  B  *  o 


-V  e-*k^x~X°^Sin0  if  m  -  o. 


Similarly 


A1  -  o 
m 


o  if  o 

-jkV  eJk(x‘x°)Sin0  if  m  -  o. 


In  the  above  expressions,  properties  of  the  Bessel  function 


are  used. 


J  (o)  o  if  m  + 

m  =  1  if  m  - 


Expressions  (B-15)  and  (B-19)  reduce  to 


V.  -  B  -  -V  Sine  eJk(x“*o^ 
in  o  o 

31  Vln  "  Tln  "  ®o  “ 


ii 


(B-20) 
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Cue  II.  Wave  Traveling  Normal  to  Axis  of  Symmetry 

If  «  =  o,  the  yx  axis  coincides  with  the  y  axis,  and  we  can  get  the  velocity  distribution  for  a  wave 
traveling  along  the  -y  axis.  If  we  want  to  consider  the  wave  traveling  along  the  +  y  axis,  we  have  to  use 
*  =  180°. 


In  the  following  we  will  consider  the  «  =  180°  case. 

A  -  -  ^  Cose  (-1)“  ejky°  J  (kr) 
m  i.  m 


A1  -  ~  k  SineCose  e  jky°  J*  (-kr) 
m  Z  m 


B  »  0 
m 

vin  •  23*i  +  Li  a"'1  (vrVi) 
vL  -  23A1  *  S-i 23  "'1 


Cosing 

Cosm$ 


(B-21) 
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Appendix  C 

Analysis  of  Scattering  of  the  Plane  Wave 
From  Rigid  and  Flexible  Structures 

The  net  component  of  the  particle  velocity  normal  to  the  surface  Vn  is  equal  to  the  sum  of  the 
normal  component  of  the  incident  velocity  and  the  scattered  velocity. 

V11  =  n  +  v®  (C-l) 

where 


Vn  =  v" n,  V^n  =  Vin  (C-2) 

Vs  =  VSn 

V11  =  vin  +  Vs. 


Scattered  components  of  the  pressure  and  the  velocity,  ps  and  Vs,  respectively,  satisfy  the  linear 
algebraic  counterpart  of  the  appropriate  boundary  value  problem2  of  the  Helmholtz  equation: 

[Hj  {Ps}  =  [G]  {V8}.  (C-3) 

Incident  velocity  and  pressure,  v'  and  p‘,  respectively,  are  related  as  follows: 


+  p*  dx 
»x 


P  9V1  -  ^P1 

f  9t  3x 


(C-4) 


For  harmonic  velocity,  equation  (C-4)  reduces  to 

Pi(x,t)  =  -cp^fx.t).  (C-5) 

The  total  acoustic  pressure  at  a  point  is  the  sum  of  the  total  scattered  pressures  and  incident 
pressure  given  by  equations  (C-3)  and  (C-5),  respectively. 


Substituting  equation  (B-23)  into  equation  (B-24),  we  get 


{p}  -  (p1)  +  [hJ_1  [g]  j  {v“}-{Vin}J 

p  --cpf{vi}+  [hJ-1  £gJ  {v11}-  ^hJ”1  £gJ  {vin}.  (C-6) 


Note  that  in  the  above  step  we  constructed  the  matrix  vectors  from  the  set  of  scalar  equations  (C-2) 
and  (C-5). 


C-l 


I 


Case  I.  Rigid  Structure 

If  the  structure  is  assumed  rigid,  the  normal  velocity  {Vn}  of  the  structure  is  zero,  and  equation 
(C-6)  reduces  to 

(p>  *  cOftV1}  -  [h]'1  [cjtV10}  .  (C-7) 

The  listing  of  the  computer  program  uses  this  set  of  equations. 

Case  II.  Flexible  Structure 

The  equations  of  motion  of  the  finite  element  model  of  the  structure  are  written  as  follows: 

[M]  {x}+[k]{X}  -  {F}+[A]  {p}.  (C-8) 

Substituting  equation  (C-6)  into  equation  (C-8)  and  noting  that  the  vector  {Vn}  is  related  to  the 
vector  { X } ,  we  can  rewrite  equation  (C-9)  as 

[M]  (x'1-[a][H]";L[g]{X}+[k]  {X}  -  {Fj-C0f  [A]  {v1>-cpf  [a][h]  _1[g]  {V111}  (C-9) 
without  introducing  the  transformation  matrix  relating  the  vectors  {X}  and  { V"}. 


TR  5821 


INITIAL  DISTRIBUTION'  LIST 


Addressee  No.  of  Copies 

CNM,  MAT-08T1  (Dr.  J.  Probus)  1 
NRL  (Dr.  Neubauer)  1 
USOC,  Code  240  1 
NAVELECSYSCOM,  PME-124  (CAPT  Cox),  ELEX  320  (Joel  Sinsky)  (2)  3 
NAVSEASYSCOM,  SEA-05H(CAPT  Jurgen)  1 
DTNSRDC  (Dr.  Wayne  Reader)  1 
DTIC  12 
DARPA  (Dt.  Theo  Kooij)  2 
SSPO  (Dr.  Phil  Selwyn)  1 


